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Abstract 


In high-altitude nuclear detonations, the proportion of pulsed X-ray energy can exceed 70%, 
making it a specific monitoring signal for such events. These pulsed X-rays can be captured using 
a Satellite-borne X-ray detector following atmospheric transmission. To quantitatively analyze the 
effects of different satellite detection altitudes, burst heights, and transmission angles on the physical 
processes of X-ray transport and energy fluence, we developed an atmospheric transmission 
algorithm for pulsed X-rays from high-altitude nuclear detonations based on scattering correction. 
The proposed method is an improvement over the traditional analytical method that only computes 
direct-transmission X-rays. The traditional analytical method exhibits a maximum relative error of 
67.79% compared with the Monte Carlo method. Our improved method reduces this error to within 
10% under the same conditions, even reaching 1% in certain scenarios. Moreover, its computation 
time is 48000 times faster than that of the Monte Carlo method. These results have important 
theoretical significance and engineering application value for designing satellite-borne nuclear 
detonation pulsed X-ray detectors, inverting nuclear detonation source terms, and assessing 


ionospheric effects. 
Keywords: High-Altitude Nuclear Detonation; Atmospheric Transmission; Pulsed X-rays; 
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1 Introduction 


Nuclear detonations can be classified into three main categories based on their location: 
underground or underwater, ground and surface, and low and high altitudes. High-altitude nuclear 
detonations differ significantly from other types of detonations [1]. According to publicly available 
data from the U.S. military [2], a typical low-altitude detonation mainly transmits energy through 
shock waves, whereas high-altitude nuclear detonations are primarily radiation driven, and the 
proportion of X-ray energy released (TNT equivalent) at the moment of detonation can reach up to 
85% of the total equivalent [3]. Pulsed X-rays formed by high-altitude nuclear detonations serve as 
specific monitoring signals for such detonations [4-6]. After detonation, the X-rays interact with the 
atmosphere and eventually reach the satellite-borne X-ray detectors. Studying the physical processes 
of atmospheric X-ray transmission is essential for monitoring high-altitude nuclear detonations and 
detecting illegal nuclear tests under the Comprehensive Nuclear Test Ban Treaty. 

Owing to treaty limitations, high-altitude nuclear weapon testing cannot be conducted. 
Therefore, current research on X-ray atmospheric transmission calculations for high-altitude nuclear 
detonations uses two primary methods: Monte Carlo and analytical methods. Harris et al. of the U.S. 
Defense Nuclear Agency developed the ATR calculation code, which uses Monte Carlo simulations 
to calculate the radiation environment generated by nuclear detonation X-rays and neutrons at 
different altitudes in the atmosphere [7]. The Ballistic Research Laboratory developed the FLAIR 
code, which obtains the photon energy, angle, and time for X-ray and gamma-ray transmissions 
using Monte Carlo database calculations [8]. Liu et al. used the Monte Carlo method to simulate the 
characteristic ring pinhole imaging of nuclear detonation X-rays under equivalent temperature 
blackbody spectra of 1.4 and 3.8 keV [9]. They established a reverse nuclear detonation equivalent 
table and roughly inferred the nuclear detonation equivalent through a table look-up method. Liu et 
al. analyzed a nuclear detonation debris cloud that generates pulsed X-rays and established a fluid 
dynamics model of debris motion from a near-space nuclear detonation with a TNT equivalent of 1 
kt-10 Mt and a burst height of 30-80 km [10]. The results showed that the maximum height, 
horizontal radius, and speed of the debris cloud increase with increasing explosion height and TNT 
equivalent. Ou-Yang et al. used analytical numerical simulation methods to investigate the effects 
of nuclear detonation pulsed X-rays on the ionization and evolution of the atmosphere at different 
burst heights and radiation angles [11]. Xu performed analytical calculations of nuclear detonation 
X-ray energy deposition and demonstrated that the energy deposition is proportional to the TNT 
equivalent of the nuclear detonation [12]. In general, Monte Carlo methods can achieve high 
computational accuracy, but they are computationally expensive and cannot be applied to scenarios 
requiring rapid responses, such as nuclear detonation detection. Furthermore, Monte Carlo methods 
cannot directly infer the TNT equivalent from the output of nuclear detonation detectors. An 
effective way for improving both the calculation speed and inversion of the nuclear detonation 
equivalent is to use analytical methods. However, the current analytical methods consider only the 


atmospheric attenuation effect on the direct X-ray transmission orientation. Owing to their failure 
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to consider scattered X-rays, analytical methods have inferior calculation accuracy than Monte 
Carlo methods. Based on Monte Carlo simulations, Xiao et al. concluded that the high-energy 
nuclear detonation X-ray fluence originates mostly from non-collision X-rays, whereas the low- 
energy fluence originates mainly from X-rays scattered during transmission [13]. High-altitude 
nuclear detonations constitute a significant proportion of low-energy X-rays. Thus, it is necessary 
to compute scattered X-rays using analytical methods. 

To address the above issue, this study established an atmospheric transmission model for 
pulsed X-rays in high-altitude nuclear detonations. The scattering effect was corrected using build- 
up factors and scattering correction coefficients. We used the proposed method to calculate the 
energy fluences of high-altitude nuclear detonation pulsed X-rays at different satellite altitudes, 


burst heights, TNT equivalents, and transmission angles. 
2 Method 


2.1 Nuclear detonation X-ray source 

The temperature of an X-ray fireball generated by a nuclear detonation can reach 10°- 107 K. 
At this temperature, the substances in the area surrounding the detonation are completely ionized, 
forming high-temperature plasma [14]. The plasma exhibits a large contrast with its surroundings 
and can be approximated as a blackbody [3, 15]. The normalized energy spectrum can be represented 


by the Planck blackbody spectrum as follows: 
15 E? 
E EOR I EON (1) 
Gr x Ty) exp(E/Ty) — 1 
where E is the energy of the X-ray photons and Ty is the equivalent temperature of the X-ray fireball. 


P(E,Ty) = 


Assuming that u = E/Ty is a dimensionless energy ratio, the normalized energy spectrum can be 
written as 
15 u’ 


PUR E oT EE 


(2) 


Some studies [16, 17] have shown that the equivalent temperature of X-ray fireballs is related 
to the design details of the nuclear weapon detonation processes and projectile materials. X-ray 
sources can be broadly divided into three categories based on the projectile type: fission bombs, 
ordinary hydrogen bombs, and enhanced X-ray weapons. The X-ray energy spectrum of a fission 
bomb can be characterized using a single blackbody spectrum with an equivalent temperature of 
approximately 1-2 keV. The X-ray energy spectra of other types of weapons can be characterized 
using composite blackbody spectra obtained by combining multiple blackbody spectra in different 
proportions. Here, the nuclear detonation X-ray source was a fission atomic bomb with an equivalent 


temperature of 1.4 keV. 


2.2 Atmospheric transmission model 
To calculate the parameters of nuclear detonation X-ray transmission to satellite-borne 
detectors, it is necessary to establish a mathematical model for nuclear detonation X-ray 


transmission through the atmosphere. The U.S. Standard Atmosphere 1976 was used in the 
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atmospheric model in this study [18]. The atmosphere above 1000 km is extremely thin and can be 
considered a vacuum environment with an atmospheric density of 0. Statistically [19], the 
atmosphere is mostly distributed below 90 km, and the ratios of major elements (C, N, O, etc.) in 
the atmosphere hardly change in the 40-90 km range. The density of the atmosphere from 90 to 120 
km is low, and the proportions of N and O vary, whereas the contents of the other components are 
very small. Monte Carlo simulations verified that the change in the proportion of each component 
in the atmosphere with altitude has a negligible effect (with a maximum relative error of less than 
0.6% compared with the results for a fixed component) on the X-ray photon energy fluence. 
Therefore, it was assumed that the atmospheric composition is independent of altitude. Supposing 
that Earth is a standard sphere, the atmospheric density is approximated to have a spherically 
symmetric distribution. 

The transmission parameters were calculated by establishing the transmission path relationship 
between detector A (where the altitude of the satellite is h4), detonation point B(B’) (where the 
burst height is hg), and Earth (where the spherical center is O, and the Earth’s radius is Rg). 6 is the 
observation angle, a is the transmission angle, and AB(AB’) is the X-ray transmission path, as 


shown in Fig. 1. 


Fig. 1. High-altitude nuclear detonation X-ray transmission 


2.3 Direct transmission 

If a nuclear weapon detonates, X-rays interact with the atmosphere surrounding Earth layer- 
by-layer, starting at the detonation point, until they reach the satellite-borne X-ray detector. Nuclear 
detonation X-rays transmitted to satellite detectors consist of two components: directly transmitted 
[20] and scattered X-rays. First, we calculated the directly transmitted X-rays. 


The atmosphere from the detonation point to the atmospheric boundary was assumed to be 
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vertically divided into N(M) layers according to the 1976 standard atmosphere. The stratification 
height of the i-th (0 <i<N(M )) atmosphere layer is AS;. The layer containing the detonation 
point is the O-th layer, with ASọ = 0. Other stratification heights 4S; are a piecewise function of 
altitude h: 

0.1,0 < h < 30 


ASi(1 <i < N(M))=}2, 30<h<90 . (3) 
5, 90 <h < 1000 


Calculation of the X-ray transmission attenuation requires the atmospheric column density ps 
along the transmission path [21]. The atmospheric column density, also known as mass absorption 
distance, is the integral of the atmospheric density per unit length over distance. Two cases were 
considered for ps according to the transmission angle a. 

I. Case 1: Transmission angle a > 90° 

This case corresponds to detonation-point location B in Fig. 1. The transmission paths are 

BC > BD > »--- > BA. The layered paths are BC > CD > -:- > A. In this case, the transmission 


paths are all above the burst height hg of point B. The total number of layers N can be written as 
N= Ni + N- j (4) 


where N; and N, represent the number of layers within the 30-90 km range and 90-1000 km range, 


respectively. 


90 — hg 


N = 30 < hg <90, (4a) 
0, else 
182, 30 < hg < 90 
—7. else 


For ACBO, according to the cosine theorem, 
OC? = BC? + OB? — 2BC x OB x cosa, (5) 


where OC = Rg + hg + AS, and OB = Rẹ + hg. Suppose that BC = X,, and 
H x cosa + |H2 x cos? q + 4(AS,? +Hx AS;) 
Xi = — Yr 6 
> (6) 
where H = 2(R; + hg), and the length AX, of the first layered path BC is 
AX, = X4- (7) 
By analogy, the recurrence relation X,,(2 < n < N) for the length of BC > BD > --. > BA is 


Hx cosa + 4 H? x cos*a+4((XiL, AS;)? +H x YL, AS;) 
x, =e ee (8) 


The length AX,,(2 < n < N) of the n-th layered path of BC > CD > --- > Ais 


AXn =Xn—Xn-1- (9) 
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The atmospheric column density of the transmission path in Case 1 can be expressed as 


N 
n-1 
Ps(hsa = 90°) = X AXnxpy(Re +) Asi), (10) 
i 
n=1 


l 


where py(h) is the atmospheric density at altitude h. 
II. Case 2: Transmission angle a < 90° 

This case corresponds to detonation-point location B’ in Fig. 1. The transmission paths are 
B'F > B'G > B'C' > B'D' --- > B'A. The layered paths are B'F > FG > GC’ > C'D'--» > A. In 
this case, some layered paths are partially above the burst height hg, such as GC’; others are below 
hg, such as B'G. To distinguish from the layers above the burst height hg, the stratification height 
in section B'G is expressed using As;, whose physical meaning and expression are the same as AS; 
in Equation (3). 

In this case, the total number of layers M can be written as 

M =N + Mp'g, (11) 

where N is given by Equation (4) in Case 1, and M g'g represents the number of layers in section 
B'G. Mp'g is related to the shortest altitude hmin of the transmission path to Earth’s ground and 


burst height hg, where 


1 
OF = (Rg + ha) x sin ô = 5 H4 X sin ô, (12) 
hmin = OF — Rp. (13) 


Layer M g'g can be divided into three parts: 
Mp'g = Mı + M, + M3, (14) 
where M1, M3, and M3 represent the number of layers below 30 km, between 30 and 90 km, and 


between 90 and 1000 km, respectively. They can be expressed as 


M, = 1 min = 30 (14a) 
0, else 
hp a 3 
5 (hmin < 30) and (30 < hg < 90) 
30, (hmin < 30) and (hg > 90) 
hg — hmi 
M, = aa eae (30 < hmin < 90) and (30 < hg < 90), (14b) 
90 — hmi 
S (30 < hmin < 90) and (hg > 90) 
0, else 
hp a 90 
=—, (hmin S 30) and (hg > 90) 
Amin — 90 
WE D (30 < hmin < 90) and (hg > 90) (140) 
hp — hmi 
ee Amin >90 
0, else 


In AB'C'O, where B'C' = X4, the length AX; of the first layered path GC’ above the burst 
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height hg is 
AX, = X,—2FB' = X, -2 (H/2)? — (H,/2)? x sin? 6. (15) 


The other recurrence relations for X, and AX,, are the same as those in Equations (8) and (9) in 
Case 1. 
For B'G below burst height hg, the detonation point can be equated to point F, where point 
B’ is equivalent to the satellite. According to the geometric relationship, 
B'G = 2B'F, (16) 
the transmission paths of section B’F have a recurrence relationship similar to that of Case 1. The 
first-layer transmission path length x, and layered path length Ax, below the burst height hg have 


the following relationship: 


Ax, =x, = y (OF + As,)? — OF? = |H; X sin ô x As, + As,”. (17) 


Similarly, the transmission path length xm and layered path length Ax,,(2 < m < Mpg'g) below 


burst height hg have a recursive relationship: 


tn = |Ha x sind x Yi, As; + (x™, Asi), (18) 


AX mi = Xm — Xai (19) 


The atmospheric column density of B'F is defined as extra ps, and is expressed as 


M 
m-1 
Ps, = > AXm X Py (hmin + A Asi) - (20) 
net i=0 


The atmospheric column density of the transmission path in Case 2 can be written as 


n-1 


N 
ps(hg,& < 90 °) = D AXn x Pv (hs + >, asi) + 2Ps, (21) 
E| i=0 


After calculating the atmospheric column density ps for Cases 1 and 2, the X-ray transmission 
attenuation by the atmosphere can be calculated based on ps. The mean free path of photons is a 
general measure of the opacity of a substance and the ability of a medium to absorb radiation [22]. 
The mean free path A of X-rays with energy E between the detector and detonation point is defined 
as 

ACE) = Um(E) X ps, (22) 
where Um(E) is the mass absorption coefficient of X-rays with energy E, which was interpolated 
based on discrete data provided on the official NIST website [23]. 

After detonation, the radiation generated by the X-ray source spreads spherically to the 

surrounding areas. The radiant intensity of nuclear detonation pulsed X-rays with energy E at the 


detector can be expressed as 


E+AE 
1(E) =| Ko g TROR UX TL) 


e 4) dE, 23 
E Amt X Ro” oy 
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where I(E) is the intensity of direct-transmission X-rays with energy E; Kg is the equivalent- 
energy conversion factor (Kg = 2.6114 x 1025 keV/kt), with the equivalent expressed in kt TNT; 
f(t) is the X-ray-normalized time spectrum; Rọ is the distance from the detonation point to the 
satellite-borne X-ray detector; and 7 is the proportion of X-ray equivalent Qx in the total detonation 
equivalent Q. ņ is related to the nuclear weapon design and its value is generally taken from 0.7 to 


0.85 [3]. The expression for the relationship is as follows: 
Qx 
y=—. (24) 

Q 

2.4 Build-up factor 

In addition to being attenuated by collisions with atmospheric particles, X-rays undergo 
scattering. In this study, the correction for direct-transmission orientation scattering was made using 
the build-up factor. Bridgman defined the build-up factor as “direct and unreacted photons plus the 
scattered contribution.” [24] The build-up of X-ray intensity J for direct and scattered transmissions 
in the direct-transmission orientation comprises three components: intensity of direct transmission 


1©), intensity after being scattered once in the direct-transmission orientation /™, and intensity 


after m(m = 2) scattering 1”. 
P= 1 4 7) +, (25) 


The single-scattering intensity was derived by Bigelow and Winfield [25], and the results of multiple 
scattering were derived by Renken using the Boltzmann transfer equation [26]. The relationship 
between the build-up factor Fg and radiation fluence in an infinite homogeneous atmosphere for X- 
rays with energy E is given by 
I(E) = 1 (E) x F(E). (26) 
Tayler fitted the build-up factor as a function [27] based on Monte Carlo simulation results: 
Fg (E) = A; X exp(c, X A(E)) + Az X exp(c2 X A(E)), (27) 
where A4, A2, C1, and c, are coefficients related to the X-ray photon energy E. 
The X-ray build-up factor is related to Compton scattering. Kalansky concluded that, based on 
the Compton scattering of photons, the value of the build-up factor can be divided into three 


segments depending on the X-ray energy E [28]: 
I. For E < 12 keV, the Compton scattering is negligible, with Fg = 1. 


IL. For 12 keV < E < 750 keV, Talyer's fitting equation can be used and the Compton 


scattering cross-section of X-rays with energy greater than 100 keV starts to decrease. 


Il. For E > 750 keV, the electron pair effect is dominant, and this fitting formula is no 
longer applicable. 


For the nuclear detonation X-ray source adopted in this study (fission atomic bomb), only 
approximately 1% of the X-ray energy was greater than 10 Ty, which corresponded to an energy of 


14 keV. The proportion of X-rays with E > 750 keV was low and could be ignored [3]. 
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The energy fluence represents the total energy of photons entering a spherical body per unit 
cross-sectional area during X-ray emission [29]. The initial energy fluence ®p is defined as the sum 
of the X-ray energy fluences transmitted and scattered in the direct-transmission orientation [30]. 


This is the integral of the X-ray intensity over time after correction for the build-up factor. 


nxQ ft E ae 
aol rae | F(E) x P(E, Tx) Xe dE. (28) 


Po = | I dt = Kg X 

When a nuclear weapon is detonated at a high altitude, the spectrum peak and duration remain 

in the order of nanoseconds to hundreds of nanoseconds [17]. Thus, it can be assumed that the X- 

ray energy spectrum does not vary with time during the observation [16]; therefore, the effect of 
time t was ignored. 

In the direct-transmission orientation, the transmission coefficient [31] TCE) represents the 

probability that an X-ray with energy E is transmitted through the atmosphere to the satellite-borne 

detector without considering geometric attenuation. Owing to the build-up factor, the sum of the 


transmission coefficients may be greater than 1. The sum of the transmission coefficients at all 


energies is expressed as 


XŒ) = È Fa(E) x PCE, Ty) x exp (-A(E)). (29) 
E=0 
By ignoring the influence of time t, the initial energy fluence can be written as 
nxQ > 
Py = Ko X ———_5X ) T(E). 30 
vaan (E) (30) 


2.5 Monte Carlo-based scattering coefficient correction 

Because the nuclear detonation source can be approximated as a point source, the X-rays 
emitted after detonation radiate spherically toward the surrounding area. Therefore, in numerical 
simulations, in addition to considering the correction of direct-transmission scattering using build- 
up factors, it is necessary to consider the scattering effect of X-rays from other orientations. This 
study established a correction factor table to correct the energy fluence based on Monte Carlo 
simulation results. X-ray transport simulations were performed using the SuperMC Monte Carlo 
simulation software. 

The average energy fluence £ of a single photon at the satellite altitude was calculated using 
Monte Carlo simulation [32]. Energy fluence ®yc can be calculated by the average energy fluence 
of individual photons £ amplifying the number of X-rays ny when the detonation equivalent is Q: 

Puc = Ny XZ, (31) 


where ny can be expressed as 


wW n x Ei 
ma > Ko x x| P(E,T, dE. (32) 


i=1 l Ei-1 


Here, W is the number of energy intervals, and the X-ray energy spectrum was divided into 


2 x 10° energy intervals. €;(1 < i < W) is the average photon energy of the i-th energy interval. 
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The detector counting results simulated using the Monte Carlo method were considered valid 
when the statistical error was within 5% [33]. The maximum permissible number of particles for 
each simulation condition was 2 X 10°. Under these limitations, the Monte Carlo simulation results 
demonstrated that the lowest credible burst height was approximately 42 km. 

By comparing the data of Monte Carlo simulated energy fluence ®yc with the data of initial 
energy fluence ®o, we found a strong linear relationship between the two. Therefore, X-ray energy 
fluence @ at satellite altitudes can be further corrected and expressed as 

$ = K(hg) X bo (33) 
where K is the scattering correction coefficient related to burst height hp (listed in Table I), and R? 
is the coefficient of determination. After testing in areas with burst heights hg > 100 km, the 
simulation results show that K ~ 1 is applicable for 100 km < hg < 1500 km. 


Table I. Scattering correction coefficients applicable to different heights 


hg / km K R? hpg / km K R? hg / km K R? 
42 1.4550 0.9999 60 .3168 0.9993 78 1.0625 0.9995 
44 1.7756 0.9999 62 .2717 0.9993 80 1.0513 0.9995 
46 1.8236 0.9999 64 .2329 0.9994 85 1.0454 0.9990 
48 1.8479 0.9999 66 .1985 0.9994 90 1.0477 0.9982 
50 1.5682 0.9999 68 .1675 0.9995 95 1.0214 0.9997 
52 1.5362 0.9999 70 .1401 0.9995 100 1.0071 0.9920 
54 1.4864 0.9996 72 1.1160 0.9995 >100 ~] - 
56 1.4275 0.9994 74 .0948 0.9995 
58 1.3690 0.9993 76 .0770 0.9995 

3 Results 


3.1 Comparison with Monte Carlo simulation results 
The calculation conditions for the analytical simulation of the atmospheric transmission of 


pulsed X-rays from high-altitude nuclear detonations are listed in Table II. 


Table II. Numerical simulation calculation conditions 


Variable Value and Unit Bullet type 
Satellite altitude h4 20000 km 
Burst height hg 42 km- 1500 km 
Transmission angle a 89°— 180° Fission atomic bomb (Ty = 1.4 keV) 
Detonation equivalent Q 100 kt TNT 
Proportion of the X-ray equivalent n 0.7 


For each burst height, we randomly selected no less than six transmission angles from 89, 90, 
95, 100, 110, 120, 130, 140, 150, and 160° and calculated the corresponding initial energy fluence 
®, corrected by the build-up factor and ® yç using the Monte Carlo method. Figure 2 shows that 
@, and ®y- exhibit a good linear relationship. The slope of the fitting line in Fig. 2 represents the 


scattering correction coefficient. 
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hg = 42~50 km (a) hg = 85~130 km (b) 
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Initial energy fluence (keV/cm?) x10" Initial energy fluence (keV/cm?)x10"° 
Fig. 2. Fitting diagram of initial energy fluence and Monte Carlo results 
The energy fluence ® was corrected by the build-up factor and scattering correction coefficient. 
For cases where the transmission angle was below 90°, the relative error between ® and Puc 
exceeded 15%, when the burst height was less than 120 km. The main reason for the large relative 
error in these cases is that the interaction process between the X-rays and atmosphere becomes more 
complex after adding the extra atmospheric column density, and the scattering process becomes 
more difficult to correct. If the transmission angle is appropriately selected based on different burst 
height conditions, the relative error is within a reasonable range (within 15% is permissible for 


detection), as listed in Table III. 


Table II. Suggested transmission angles for different burst heights 


hg / km Suggested availability range of æ 
42~64 2> 120° 
64~74 => 110° 
74~90 => 100° 
90~100 295° 
100~130 2> 90° 
Above 130 => 89° 


To validate the algorithm, we compared the proposed analytical method and another method 
without scattering correction with the Monte Carlo method. The results are presented in Table IV, 
where the uncorrected energy fluence ¢ is the energy fluence at the detector for direct transmission 


of attenuated X-rays. 
Table IV. Energy fluence calculated using the Monte Carlo method, uncorrected analytical method, and proposed method 


Energy fluence (ke V/cm’) 


Monte Carlo Uncorrected Relative error Initial energy Relative error Relative error 
hg / km ale Energy fluence 
energy fluence energy fluence between @ and fluence between ®y © between ® and 
Puc p Duc Do and mc Puc 

120 2.3785x10° 7.6610x10° -67.79% 1.5632x10° -34.28% 2.2744x10° -4.38% 

42 130 6.2691x10° 2.3916x10° -61.85% 4.3270x10° -30.98% 6.2959x10° 0.43% 

140 1.2065x10" 4.9709x10° -58.80% 8.3617x10° -30.69% 1.2166x10’ 0.84% 

160 2.5130x10" 1.100410’ -56.21% 1.7199x10" -31.56% 2.502410" -0.42% 

110 7.8992x10" 3.7116x10’ -53.01% 5.0120x10" -36.55% 7.4498x10" -5.69% 

54 120 2.0101x108 1.0322x108 -48.65% 1.2743x108 -36.60% 1.8942x108 -5.77% 

140 5.3173x108 3.1004x108 -41.69% 3.5159x108 -33.88% 5.2260x108 -1.72% 

150 6.9533x108 4.2153x108 -39.38% 4.6872x10° -32.59% 6.9670x108 0.20% 

78 100 5.1513x10° 4.3856x10° -14.86% 4.4090x10° -14.41% 4.6845x10? -9.06% 
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120 1.2748x10'° 1.1749x10'° -7.84% 1.1761x10'° -7.14% 1.2496x10'° -1.98% 
130 1.5688x10'° 1.4671x10'° -6.48% 1.4682x10'° -6.41% 1.5599x10'° -0.57% 
140 1.8193x10'° 1.7172x10"° -5.61% 1.7182x10'° -5.56% 1.8256x10'° 0.35% 
90 1.4196x10"° 1.3528x10'° -4.71% 1.3533x10" -4.68% 1.3629x10!° -4.00% 
100 95 2.2652x10"° 2.1767x10"" -3.91% 2.1769x10'° -3.90% 2.1924x10'° -3.22% 
100 2.6434x10" 2.5813x10"° -2.35% 2.5814x10" -2.35% 2.5997x10!° -1.65% 
110 3.1234x10"° 3.0919x10!° -1.01% 3.0920x10'° -1.01% 3.1139x10"° -0.30% 
Based on the data in Table IV, the results obtained by the uncorrected analytical method are 
significantly different from those obtained by the Monte Carlo simulation. Because the uncorrected 
analytical method does not consider the scattering effect, all its energy fluence values are lower than 
those of the Monte Carlo simulation. Under a burst height of 42 km and transmission angle of 120°, 
the relative error between the two reaches 67.79%. An improvement in accuracy is evident after 
scattering correction, with an increase of more than 60% in this case. 
To further verify the accuracy, the results under different satellite altitudes, burst heights, and 
equivalent conditions were compared, as listed in Table V. 
Table V. Verification of randomly selected different detonation equivalents and altitudes 
Energy fluence (ke V/cm’) 
Relative Relative Reaive cier 
Q/ktTNT h,/km hg/km al’? error error 
Puc Q Po p between ® 
between p between ®y 
and ®yc 
and mc and ®yc 
120 2.9410x10° 1.9120x10° -34.99% 2.0616x10° -29.90% 2.7147x10° -7.69% 
10000 60 130 4.8608x10° 3.3776x10° -30.51% 3.5544x10° -26.88% 4.6805x10° -3.71% 
150 8.8068x10° 6.5290x 10° -25.86% 6.7381x10° -23.49% 8.8727x10° 0.75% 
160 1.0376x10'° 7.8116x10° -24.71% 8.0298x10° -22.61% 1.0574x10'° 1.91% 
120 5.4753x10° 3.8447x10° -29.78% 4.0021x10° -26.91% 5.0118x10° -8.46% 
10000 63 130 8.5729x10° 6.3309x10° -26.15% 6.5018x10° -24.16% 8.1422x10° -5.02% 
140 1.1731x10"° 8.9292x10° -23.88% 9.1088x10° -22.35% 1.1407x10'° -2.76% 
100 150 1.4639x10'° 1.1358x10'° -22.42% 1.1544x10'° -21.15% 1.4456x10'° -1.25% 
120 2.0115x10" 6.2783x10° -68.79% 1.0734x10" -46.64% 1.9317x10" -3.97% 
15000 45 130 4.7660x10" 1.6490x10’ -65.40% 2.5733x107 -46.01% 4.6309x10’ -2.84% 
150 1.2675x108 4.703110’ -62.89% 6.7138x107 -47.03% 1.2082x108 -4.67% 
160 1.6430x108 6.1871x10" -62.34% 8.6410x107 -47.41% 1.5550x108 -5.35% 
110 4.0104x10° 1.2709x10° -68.31% 2.1608x10° -46.12% 3.9930x10° -0.43% 
36000 48 120 1.2648x10" 4.7400x10° -62.52% 7.0798x10° -44.02% 1.3083x10’ 3.44% 
140 4.1604x10’ 1.6980x10’ -59.19% 2.262210" -45.63% 4.1802x10’ 0.48% 
160 7.0581x10" 3.0006x10" -57.49% 3.8058x10" -46.08% 7.0328x10" -0.36% 
110 1.4151x108 6.2294x10" -55.98% 8.7038x107 -38.49% 1.3154x108 -7.05% 
18000 53 120 3.7582x108 1.8078x108 -51.90% 2.3010x108 -38.77% 3.4776x108 -7.46% 
130 6.9098x108 3.5612x108 -48.46% 4.2852x108 -37.98% 6.4762x108 -6.28% 
200 140 1.0438x10° 5.6740x108 -45.64% 6.5904x10° -36.86% 9.9600x108 -4.58% 
110 1.5772x10"° 1.2566x10'° -20.32% 1.2706x10'° -19.44% 1.4333x10'° -9.13% 
13000 71 120 2.5274x10" 2.1047x10!° -16.72% 2.1171x10" -16.23% 2.3882x10'° -5.51% 
130 3.4483x10'° 2.9472x10"° -14.53% 2.9587x10"° -14.20% 3.3376x10'° -3.21% 
140 4.2982x10!° 3.7358x10"° - 13.08% 3.7467x10'° -12.83% 4.2265x10!° -1.67% 
100 1.7172x10" 1.5166x10'° -11.69% 1.5188x10'° -11.55% 1.5914x10" -7.33% 
110 2.5966x10"° 2.3960x10!° -7.13% 2.3974x10" -7.67% 2.5119x10" -3.27% 
500 36000 83 
120 3.2144x10" 3.0240x10!° -5.92% 3.0250x10'° -5.89% 3.1695x10'° -1.40% 
160 4.5756x10'° 4.4064x10'° -3.70% 4.4070x10'° -3.69% 4.6175x10'° 0.91% 


The data in Table V indicate that the relative error between ® and ®yc can be controlled 


within 10% by interpolating the scattering correction coefficient K for burst heights that are not 


listed in Table I. Moreover, the K value applies to satellite altitudes greater than 10000 km. Good 


agreement between the proposed analytical method and Monte Carlo simulation can be observed 
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for the selected available transmission angles. 

The computation time [34] of the proposed analytical algorithm is reduced to 1/48000 of the 
time taken by the Monte Carlo method. Table VI presents a comparison of the time required by the 
two methods to calculate six different burst heights, with each height corresponding to 10 


transmission angles in parallel. 


Table VI. Comparison of calculation time between the Monte Carlo and proposed method 


Calculation ; Number of Calculation Calculation . 
simulated photons at . Calculation platform 
method ; method time 
each burst height 
c i System: Windows7 
oncurren 
Monte 9 y CPU: Intel Xeon Gold 6254 (single 18 cores and 36 
2x10 (MPI 3610 minutes 
Carlo thod) threads, 2pcs, total 36 cores and 72 threads) 
— RAM: 512GB 
System: Windows 10 
Analytical - Sequential 4.5 seconds CPU: Intel Core i17-8650U (4 cores and 8 threads) 
RAM: 8GB 
3.2 Numerical simulation calculation results 
— 50km 
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Fig. 3. Variation of atmospheric column density 

First, the variation in atmospheric column density corresponding to different transmission 
angles for different burst heights was calculated, as depicted in Fig. 3. As the transmission angle 
decreases, the transmission path and atmospheric column density increase. The atmospheric column 
density changes considerably when the transmission angle is approximately 90°. This is because 
when the transmission angle is less than 90°, the X-rays pass through the low-altitude and high- 
density atmosphere below the detonation point, resulting in a significant increase (2 pgs, is added) 
in the atmospheric column density. At different burst heights, even at the same transmission angle, 
the atmospheric column density differs by one order of magnitude. However, for transmission 
angles greater than 100°, the decrease in atmospheric column density exhibits a similar pattern of 


variation. 
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Fig. 4. Build-up factors for different burst heights and transmission angles 

Subsequently, the X-ray build-up factors were calculated for different burst heights and 
transmission angles, as displayed in Fig. 4. The build-up factor increases from 0 to 75 keV before 
gradually decreasing. The burst height and transmission angle significantly affect the build-up factor. 
The build-up factor varies substantially for transmission angles above and below 90°. The 
atmospheric density is relatively high when the burst height is less than 80 km and the build-up 
factor is relatively large. It can be concluded that there are many scattered X-rays in the direct- 
transmission orientation. Moreover, X-rays in the energy band of 50-150 keV have scattering 
intensities exceeding direct transmission by a factor of 100 in the case of a 50 km burst height and 
89° transmission angle. As the burst height increases, the atmosphere gradually becomes thinner 
and the scattering effect decreases. The build-up factor is close to 1 above a burst height of 100 km, 
and the influence of direct-transmission orientation scattering can be ignored. 

Finally, as illustrated in Fig. 5, under an initial explosive equivalent Q = 100 kt TNT and X- 
ray equivalent share of 0.7, the burst height of 80 km is considered as the variation boundary. This 
is because the X-ray mean free path is approximately equal to the atmospheric elevation at that 
altitude [16]. Therefore, for detonations above 80 km, X-rays emitted upward escape from the 
atmosphere, whereas X-rays with burst heights below 80 km are strongly absorbed during 
transmission, and the entire energy emitted downward is deposited in the atmosphere. The energy 


fluence spans four orders of magnitude. At a burst height of 42 km and transmission angle of 120°, 
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the energy fluence is 2.2744 x 10° keV/cm?. However, it can reach 5.9504 x 10!° keV/cm’ at a burst 
height of 1500 km and transmission angle of 180°. These results provide a solid basis for the design 


of satellite-borne X-ray detectors. 
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Fig. 5. Energy fluence for different burst heights and transmission angles 
4 Conclusion 


In this study, we developed an atmospheric transmission algorithm for pulsed X-rays from 
high-altitude nuclear detonations based on scattering correction. This method performs scattering 
correction on the energy fluence of nuclear-detonation pulsed X-rays transmitted through the 
atmosphere in the altitude range of 42-1500 km to a satellite at 10000-36000 km (altitude for 
geosynchronous satellites). We provided the recommended applicable range of the transmission 
angle, and the data showed good agreement between the proposed analytical method and Monte 
Carlo simulation for the selected available transmission angles, which can meet the requirements of 
high-altitude nuclear detonation monitoring. The proposed method reduced the time cost to 1/48000 
of that of the Monte Carlo method. The maximum relative error between the simulation results of 
the traditional analytical method and Monte Carlo method was 67.79%. Using the proposed method, 
this error could be controlled to within 10% under the same calculation conditions, and even within 
1% under certain conditions. The value of the scattering correction coefficient in the indirect- 
transmission orientation decreased gradually with an increase in burst height. Scattering was 
concentrated in the region below a burst height of 100 km, and a burst height of 80 km was 
considered the variation boundary of the energy fluence. The proposed method has great theoretical 
significance and engineering application value for the design of satellite-borne X-ray detectors, 


inversion of the nuclear detonation equivalent, and assessment of ionospheric effects. 
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